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ABSTRACT 

The  delivery  truck  problem  is  one  in  which  a  truck  is 
loaded  with  m  packages,  one  package  to  be  delivered  to  each 
of  m  destinations.   The  amount  of  fuel  consumed  by  the  truck 
is  directly  dependent  upon  the  current  total  weight  of  the 
truck,  which  includes  both  the  weight  of  the  packages  and 
the  amount  of  fuel  remaining  in  the  tank.   The  problem  is 
to  determine  a  sequence  in  which  to  deliver  all  m  packages 
which  will  minimize  total  fuel  consumption.   A  branch  and 
bound  algorithm  for  obtaining  optimal  solutions  to  the 
delivery  truck  problem  is  presented,  along  with  several 
sample  problems  with  their  solutions.   A  brief  report 
of  computational  experience  is  included. 
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I.   INTRODUCTION 

The  delivery  truck  problem  is  one  in  which  a  truck  is 
loaded  with  m  packages,  one  package  to  be  delivered  to  each 
of  m  destinations.   The  amount  of  fuel  consumed  by  the 
truck  is  directly  dependent  upon  the  current  total  weight 
of  the  truck,  which  includes  both  the  weight  of  the  packages 
and  the  amount  of  fuel  remaining  in  the  tank.   The  problem 
is  to  determine  a  sequence  in  which  to  deliver  all  m  packages 
which  will  minimize  total  fuel  consumption. 

This  problem  is  closely  related  to  the  traveling 
salesman  problem.   It  differs  in  that  the  route  is  not 
closed  and  that  fuel  consumed  on  each  leg  of  the  journey 
is  not  only  dependent  upon  the  length  of  the  leg,  but  also 
the  current  weight  of  the  truck. 

This  thesis  is  the  presentation  of  an  algorithm  for 
obtaining  optimal  solutions  to  the  delivery  truck  problem. 
The  methodology  of  branch  and  bound  is  used  to  implicitly 
enumerate  all  possible  solutions.   Several  example  problems 
with  solutions  are  presented  followed  by  a  brief  report 
of  computational  experience  for  problems  with  m  ranging 
from  five  to  twenty. 


II.   RELATED  RESEARCH 

The  delivery  truck  problem  is  clearly  a  problem  in 
the  area  of  discrete  mathematical  optimization.   Reference 
8  is  a  current  summary  of  results  applicable  to  such 
problems . 

Current  research  on  the  ordinary  traveling  salesman 
problem  is  summarized  in  Ref.  2,  in  which,  after  an  analysis 
of  computational  results,  it  is  recommended  that  dynamic 
programming  [Ref.  5]  be  used  for  problems  with  thirteen 
cities  or  less,  and  that  branch  and  bound  [Ref.  9]  be  used 
for  the  rest. 

An  algorithm  for  obtaining  optimal  solutions  to  the 
modified  traveling  salesman  problem  in  which  no  return  to 
the  starting  city  is  required  is  presented  in  Ref.  7. 

In  an  extension  of  the  traveling  salesman  problem  called 
"The  Delivery  Problem"  [Ref.  4]  more  than  one  salesman  is 
available  to  visit  the  predetermined  set  of  cities. 

Another  extension  of  the  traveling  salesman  problem, 
sometimes  referred  to  as  the  "Milk  Bottle  Problem",  [Ref.  4] 
is  that  in  which  the  milkman,  (salesman),  is  to  deliver 
one  or  more  bottles  of  milk  to  a  known  set  of  houses, 
(cities).   He  must  walk  and  wishes  to  minimize  the  amount 
of  work  expended  carrying  the  milk  bottle's.   The  problem 
is  very  closely  related  to  the  delivery  truck  problem, 
the  difference  being  that  no  account  is  taken  for  fuel 


used  along  the  way.   A  branch  and  bound  algorithm  for  the 
milk  bottle  problem  is  presented  in  Ref.  3.   This  algorithm 
is  a  direct  extension  of  the  work  of  Little,  et  al . ,  [Ref. 
6]. 

The  delivery  truck  problem  can  be  formulated  as  an 
investigation  problem.   Current  results  in  Investigation 
Theory  are  contained  in  Ref.  1.   The  algorithm  presented 
in  this  thesis  uses  an  approach  similar  to  that  of  the 
branch  and  bound  algorithm  presented  in  Ref.  1. 


III.   MATHEMATICAL  DESCRIPTION 


The  following  notation  will  be  used.   Let 


(d. . )      =  a  matrix  of  distances  between  destination 

i  and  destination  j,  i  =  l,2,...,n, 

3    =  2,3,. ..,n 
tt  =  (it,,  Tip,...,  it  )  be  a  permutation  of 

the  integers  l,2,...,n  representing  an 

order  in  which  destinations  l,2,...,n 

are  visited, 
f  =  the  amount  of  fuel  required  if  path  it 

is  followed, 
f  =  the  amount  of  fuel  required  if  destinations 

tt,  ,  iTp,  tt-,,...,  tt.  are  visited  in  that 

order.   Clearly  f  =  f   . 

TT 

F   (i,f   )  =  fuel  required  to  travel  from  destination 

IT,    '  IT  4 

3         i 

tt  .  to  destination  tt  .  given  that  i  destina- 
i  3 


tions  have  already  been  visited  and  f 
fuel  has  been  consumed. 


*i 


Let  f   =0.   Then 

TT 

o 


i    J=0      j+1      j 


The  delivery  truck  problem  is  to  select  that  path  tt  which 
minimizes  f  . 


The  algorithm  contained  herein  assumes  it  is  known 
which  destination  will  be  the  first  visited,  and  the  problem 
starts  with  the  truck  initially  located  there.   If  the 
first  destination  is  not  known,  the  truck  can  be  assumed 
to  be  located  at  a  dummy  location  added  to  the  given  set 
of  destinations,  with  travel  allowed  to  any  valid  destination 
at  zero  fuel  consumption. 

Given  that  the  truck  is  initially  located  at  destination 
one,  there  are  (n-1)!  possible  solutions. 

In  different  variations  of  the  delivery  truck  problem, 
different  fuel  consumption  functions  F  are  appropriate.  In 
this  thesis,  P  is  taken  to  have  the  following  form: 


F   (i,  f   )  -  c, (c0  -  lc£  -  f   )d 

7T.    '    TT.        12        3      7T.7T.7r. 
J         1  i     1  j 


where 


c, .=  an  arbitrarily  chosen  constant 

Cp  =  the  fully  loaded  weight  of  the  true 


k 


:_  =  the  weight  of  package  k,   k  =  l,...,n 


Note  that  if  f   =0  for  all  i,  the  problem  reduces 
7T  .  '      ^ 

to  the  milk  bottle  problem  of  Ref.  3-   If,  in  addition, 
c_  =  0,   i  =  l,...,n,  the  problem  reduces  to  the  modified 
traveling  salesman  problem  in  which  no  return  to  the 
starting  city  is  required. 
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Figure  3.1 


For  the  experimental  testing  discussed  in  Section  IV(E) 
the  constants  chosen  were  c,  =  .  004>  c?  =  8000  and 
c„  =  270,  for  i  =  i,2,...,n.   The  resulting  function  is 
used  to  illustrate  variation  in  fuel  consumption  due  to 
truck  loading  as  shown  in  Figure  3.1«   Line  a  represents 
consumption  with  i  =  0  and  f  =  0,  line  b  for  i  =  5,  f  =  0, 
and  line  c  for  i  =  5  and  f  =  500. 


IV.   THE  ALGORITHM 

This  section  begins  by  describing  how  the  branch  and 
bound  method  is  applied  to  the  delivery  truck  problem, 
followed  by  a  general  statement  of  the  algorithm.   (A  flow 
diagram  of  the  algorithm  is  contained  in  Appendix  A  and 
the  computer  program  immediately  afterwards.)   Several 
sample  problems  are  used  to  illustrate  the  output  of  the 
algorithm.   A  brief  summary  of  computational  experience 
is  presented. 

A.   APPLICATION  OP  BRANCH  AND  BOUND  TO  THE  DELIVERY 
TRUCK  PROBLEM 

The  branch  and  bound  method  consists  of  selectively 
partitioning  the  set  of  all  feasible  solutions  and  computing 
bounds  on  the  objective  function  for  each  element  of  the 
partition.   This  process  is  continued  until  one  element 
of  the  partition  containing  a  single  solution  is  obtained 
for  which  the  associated  bound  is  at  least  as  good  as  that 
of  any  other  element. 

A  branching  tree  is  initiated  and  extended  such  that 
each  node  at  the  end  of  a  branch  of  the  tree  corresponds 
to  a  set  of  solutions.   The  set  of  all  such  nodes  specifies 
a  partition  of  the  solution  space.   The  bounds  associated 
with  each  element  of  the  partition,  and  therefore  with 
each  node  at  the  end  of  a  branch  of  the  tree,  represent  a 


10 


lower  bound  on  the  amount  of  fuel  used  if  any  solution 
contained  in  that  element  of  the  partition  is  followed. 

Associated  with  every  node  of  the  tree  is  a  destination. 
The  initial  node,  representing  all  solutions,  corresponds 
to  destination  number  one,  the  initial  location  of  the 
truck.   Branching  corresponds  to  selecting  a  node  from 
which  to  branch  and  in  specifying  a  destination  for  use 
in  extending  the  path  specified  by  that  branch.   The  branch 
corresponding  to  node  one  naturally  is  selected  initially 
for  extension.   This  branch  is  extended  to  the  right  and 
the  left  by  adding  two  new  nodes  to  the  tree.   See  Figure 
4.1.   The  node  labeled  k  represents  all  solutions  in 


Figure  4.1   Initial  branching 

which  the  truck  starts  at  destination  one  and  proceeds  first 
to  destination  k.   The  node  labeled  -k  represents  all  solutions 
in  which  the  truck  does  not  go  to  destination  k  first. 
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Associated  with  each  node  at  the  end  of  a  branch  of 
the  tree  is  a  path,  starting  with  destination  one,  and 
including  each  destination  number  corresponding  to  a 
positively  labeled  node  in  the  order  in  which  they  appear 
in  the  branch. 

When  a  branch  of  the  tree  is  selected  for  extension, 
the  path  specified  by  that  branch  is  called  the  current 
path. 

The  bound  on  each  ending  node  with  a  positive  label 
applies  to  all  solutions  starting  with  the  path  specified 
by  that  node.   The  bound  on  those  with  negative  labels 
corresponds  to  all  solutions  starting  with  the  path  specified 
which  do  not  include  as  their  next  element,  all  negative 
labels  of  nodes  at  the  end  of  that  branch.   In  illustration 
of  this,  consider  Figure  *J .  2 .   The  bound  on  node  labeled  5 


Figure  k.2      Sample  tree 
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applies  to  all  solutions  starting  with  elements  1,  35  5,  in 
that  order.   The  bound  on  the  node  labeled  -5  applies  to 
all  solutions  starting  with  elements  1,  3,  and  not  having 
7  or  5  as  the  next  element. 

After  a  branch  is  selected  for  extension,  it  is  necessary 
to  determine  which  destinations  are  eligible  for  use  in 
extending  the  current  path.   This  is  done  by  determining 
all  destinations  not  in  the  current  path  and  not  prohibited 
due  to  negatively  labeled  nodes  at  the  end  of  the  branch. 

The  order  of  an  element  of  the  partition  reduces  to 
one  when  there  are  no  destinations  eligible  for  use  in 
extending  the  branch  corresponding  to  that  element.   When 
such  an  element  is  located  with  a  bound  less  than  or  equal 
to  the  bound  on  all  others,  the  path  corresponding  to  the 
sole  solution  contained  in  that  element  is  optimal. 

B.   GENERAL  STATEMENT  OF  THE  ALGORITHM 

Specific  methods  of  branching  and  bounding  are  ignored 
temporarily  while  the  general  structure  of  the  algorithm 
is  presented. 

Algorithm: 

1.  Create  the  initial  node  corresponding  to  destination 
number  one.   Compute  a  bound  for  this  node.   Go  to  2. 

2.  Select  a  branch  to  extend.   This  branch  determines  the 
current  path.   Go  to  3. 

3.  Determine  the  set  E  of  destinations  eligible  for  use 

in  extending  the  current  path.   Put  into  E  all  destinations 
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not  in  the  current  path.   Remove  from  E  all  indices  corres- 
ponding to  negatively  labeled  nodes  at  the  end  of  the  branch 
being  extended.   If  |E|  =  0,  stop.   If  not,  go  to  4. 

4.  Select  from  E  a  destination  k  for  use  in  extending  the 
current  path.   Extend  the  branch  to  the  left  and  right 
creating  two  new  nodes  labeled  k  and  -k.   Go  to  5. 

5.  Compute  bounds  for  branches  ending  with  nodes  labeled 
k  and  -k.   Go  to  2. 

Upon  termination  the  current  path  is  optimal. 

C.   BRANCHING  AND  BOUNDING  RULES 
1.   Branching 

Branching  is  a  two  part  operation,  consisting  first 
of  selecting  a  branch  to  extend,  and  next  of  selecting  a 
destination  for  use  in  extending  that  branch.   It  has 
proved  to  be  advantageous  in  branch  and  bound  to  use  as 
the  first  of  these  steps  the  selection  of  that  branch 
whose  associated  bound  is  smallest  and  specifying  tie 
breaking  rules  in  the  event  they  occur.   This  rule  is  used 
in  the  algorithm. 

The  second  part  of  the  branching  operation  selects 
a  destination  for  use  in  extending  the  current  path.   The 
algorithm  selects  that  eligible  destination  which  is  closest 
to  the  last  element  in  the  current  path,  or,  in  other  words, 
the  unvisited  destination  closest  to  the  truck's  current 
position. 
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2 .   Bounding 

The  bound  associated  with  each  branch  specifies  a 
lower  bound  on  the  amount  of  fuel  used  if  the  path  specified 
by  that  branch  is  followed.   In  illustration  of  the  method 
used,  suppose  branching  is  being  carried  out  from  a  node 
labeled  j  to  destination  k,  and  that  there  are  r  destinations 
in  the  current  path.   The  exact  amount  of  fuel  required  to 

follow  the  current  path  is  f . ,  and  the  fuel  required  to 

%  J 

travel  from  destination  j  to  destination  k  is  F.  (r,  f.). 

k     J 

The  lower  bound  on  the  node  labeled  k  is 


f  *   Fk(r,  f  )  +L 


where  L  is  a  lower  bound  on  the  fuel  required  to  visit  the 
remaining  n-r-1  destinations. 

L  is  obtained  by  considering  the  upper  triangle  of 
the  reduced  matrix  (d..)  where  i  and  j  are  not  destinations 
in  the  current  path,  selecting  the  smallest  element  in 
row  k  and  the  remaining  n-r-2  smallest  elements,  ordering 
these  distances  in  non-decreasing  order,  and  computing  the 
amount  of  fuel  required  to  traverse  these  distances  in 
the  prescribed  sequence. 

Only  the  upper  triangle  need  be  considered  since 
(d  . )  Is  symmetric  and  travel  between  destination  i  and 
destination  j  is  carried  out  at  most  once. 

It  can  be  seen  that  L  Is  a  valid  lower  bound  by 
noting  that  after  leaving  destination  k,  the  truck  must 
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traverse  an  additional  n-r-1  legs  in  its  required  journey. 
Clearly  the  shortest  distance  it  will  be  required  to  haul 
n-r-1  packages  is  the  smallest  element  in  row  k  of  the 
reduced  matrix.   A  lower  bound  on  the  next  leg  is  clearly 
the  fuel  required  to  haul  n-r-2  packages  the  smallest 
distance  in  the  reduced  matrix,  excluding  the  element  used 
for  the  first  leg.   This  procedure  is  continued  until  all 
n-1  legs  of  the  journey  are  accounted  for. 

When  the  bound  is  computed  for  the  node  labeled  -k, 
the  procedure  is  the  same  except  that  selection  of  a  distance 
for  the  first  additional  leg  after  leaving  destination  j 
is  not  restricted  to  come  from  row  k  of  the  reduced  matrix. 

D.   PROOF  OP  OPTIMALITY 

Let  S  be  the  set  of  all  solutions  to  the  problem,  i.e., 

S  =  {  tt  |  tt  is  a  permutation  of  integers  l,2,...,n}  , 

and  let  P  be  a  partition  of  S  such  that 


r      Is   2  s  ' *  *  »   m   * 


Let  b.  be  a  lower  bound  on  all  solutions  contained  in 
1 


r . ,  i.e., 


b.  <  f   for  all  tt  in  P.  . 
1  —  it  i 
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In  Step  2  of  the  algorithm  a  branch  of  the  tree  is 

selected  for  extension.   Selecting  a  branch  for  extension 

corresponds  to  selecting  an  element  of  the  current  partition 

for  refinement.   If  element  P.  is  selected  and  |P.|  =  1, 

no  further  refinement  is  possible,  and  in  Step  3>  |e|  =  0. 

At  this  point  all  n  destinations  are  contained  in  the 

current  path,  and  bound  b.  is  no  longer  a  lower  bound,  but 

j 

exact.   Note  that  the  method  of  selection  used  for  P.,  as 

J 

described  in  Section  IV(C)1,  is  to  select  partition  element 
j  such  that 


b.  =  min  b.   for  all  P,  in  P. 
j        k  k 


Hence,  P.  has  an  exact  bound  which  is  less  than  or  equal 
to  the  lower  bounds  associated  with  all  elements  of  P. 
The  single  solution  it  in  P .  is  thus  a  solution  which 
minimizes  f  . 

7T 

E.   SAMPLE  PROBLEMS  AND  SOLUTIONS 

The  following  are  sample  problems  with  optimal  solutions 
The  x  and  y  coordinates  of  the  destinations  were  selected 
from  a  uniform  distribution  over  the  interval  (0,  10). 
Coordinates  were  adjusted  slightly  when  necessary  to 
insure  that  each  pair  of  destinations  was  at  least  a  unit 
distance  apart.   The  truck  is  assumed  to  be  initially 
located  at  destination  one.   The  path  shown  is  the  one  which 
visits  each  of  the  remaining  destinations  and  uses  the 
minimum  amount  of  fuel. 
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Sample  Problem    1 
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Sample  Problem   2 
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Sample  Problem   3 


20 


F.   COMPUTATIONAL  EXPERIENCE 

Computational  experience  with  the  algorithm  is  limited, 
but  an  investigation  was  made  into  the  way  the  time  required 
to  solve  an  n  destination  problem  varies  with  n.   Seventeen 
sample  problems  were  developed  as  described  in  the  preceeding 
section;  five  with  n  equal  to  five,  six  with  n  equal  to 
ten,  and  six  with  n  equal  to  fifteen.   The  average  times 
required  to  solve  each  set  of  problems,  along  with  the 
average  number  of  bounds  computed,  is  shown  in  Table  4.1. 


n 

average 

average  # 

time  (sec . ) 

of  bounds 

5 

.38 

5.1 

10 

1.6 

145.0 

15 

99.7 

2525.0 

Table  4.1 

Although  these  results  are  sparse,  this  data  suggests 
that  time  required  increases  exponentially  with  n,  as  seen 
in  Figure  4.9.   The  circles  show  the  data  points  of  Table 
4.1.   The  dots  show  times  required  to  solve  sample  problems 
having  n  equal  to  20. 
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APPENDIX  A 
FLOW  DIAGRAM  OF  COMPUTER  PROGRAM 

The  following  notation,  not  previously  defined,  is  used 

LN       =  the  last  node  created 

HN(I)    =  destination  number  associated  with  node  I, 

for  I  =  1,  2,  . . . ,  LN 
PRED(I)   =  the  node  immediately  node  I  in  the  branching 

tree,  for  I  =  2,  . . . ,  LN 
BFN      =  node  selected  for  extension 
BTH      =  destination  selected  for  use  in  extending 

the  current  path 
BB(I)     =  bound  on  node  I 
B(I)  =  BB(I)  is  node  I  is  at  the  end  of  a  branch  in  the 

tree,  and  equals  M  otherwise,  where  M  is 

large. 
NPL      =  number  of  packages  currently  delivered 
MN       =  destination  at  which  the  truck  is  currently 

located 
F(X,Y,Z)  =  fuel  required  to  travel  distance  Z  if  X 

packages  have  currently  been  delivered  and 

Y  fuel  consumed  doing  so. 
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